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Abstract 

The work fluctuations of an oscillator in contact with a heat reservoir and driven out 
of equilibrium by an external force are studied experimentally. The oscillator dynamics is 
modeled by a Langevin equation. We find both experimentally and theoretically that, if 
the driving force does not change the equilibrium properties of the thermal fluctuations of 
this mechanical system, the free energy difference AF between two equilibrium states can be 
exactly computed using the Jarzynski equality (JE) and the Crooks relation (CR) DIUE], 
independently of the time scale and amplitude of the driving force. The applicability limits for 
the JE and CR at very large driving forces are discussed. Finally, when the work fluctuations 
are Gaussian, we propose an alternative method to compute AF which can be safely applied, 
even in cases where the JE and CR might not hold. The results of this paper are useful to 
compute AF in complex systems such as the biological ones. 

1 Introduction 

A precise characterization of the dynamics of mesoscopic and macroscopic systems is a very impor- 
tant problem for applications in nanotechnologies and biophysics. While fluctuations always play a 
negligible role in large systems, their influence may become extremely important in small systems 
driven out of equilibrium. Such is the case for nanoengines and biological processes, where the 
characteristic amount of energy transferred by fluctuations can be of the same order as that which 
operates the device. In these small out of equilibrium systems, dominated by thermal fluctuations, 
a precise estimation of the free energy difference AF between two equilibrium states A and B is 
extremely useful to increase our knowledge of the underlying physical processes which control their 
dynamical behaviour. It is well known that AF can be estimated by perturbing a system with an 
external parameter A and by measuring the work W done to drive the system from A to B. When 
thermal fluctuations cannot be neglected W is a strongly fluctuating quantity and AF may in 
principle be computed using the Jarzynski equality (JE) and the Crooks relation (CR) j2] (see 
section 2 for their definition) . Indeed the JE and CR take advantage of these work fluctuations 
and relate the AF to the probability distribution function (pdf) of the work performed on the 
system to drive it from A to B along any path 7 ( either reversible or irreversible) in the system 
parameter space. Numerous derivations of the JE and CR have been produced El El HOI 
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and JE and CR have been used in order to estimate the free energy of a stretched DNA molecule 
Q] • However the derivation of JE in the case of strongly irreversible systems has been criticized 
in ref.pl]. For this reason in a recent letter ^2], we have experimentally checked JE and CR on a 
very simple and controlled device in order to safely use these two relations in more complex cases 
as the biological and chemical ones, where it is much more difficult (almost impossible) to verify 
the results with other methods. Our system is a mechanical oscillator, in contact with a heat bath, 
which is driven out of equilibrium, between two equilibrium states A and B, by a small external 
force. In ref.^Sj we have shown that the JE and CR are experimentally accessible and valid for 
all the values of the different control parameters. However one may wonder whether the driving 
force's switching rate and its amplitude were fast and large enough. Therefore in this paper we 
extend the measurements of ref.|12j till the limits of the experimental set up. We find that even 
in these extreme cases the AF is correctly estimated by JE and CR. We also observe that these 
experimental results can be exactly derived from a Langevin equation if one takes into account 
that, experimentally, the properties of the thermal noise are not changed by the external driving 
force. These conditions are close to those pointed out in Ref. pTj for the validity of the JE, thus 
they do not fully alight the theoretical debate. However the fact of having driven the system at 
large driving amplitudes and the study of the corresponding Langevin dynamics allow us to clearly 
understand the experimental limits of applicability of the JE and CR. Independently on the theo- 
retical problems that the derivation of JE may rise, these limits are of statistical nature and they 
strongly reduce the possibility of using JE and CR when the driving energy is much larger than 
the thermal energy. As a consequence these intrinsic statistical limits make the experimental test 
of the criticisms arisen in Ref. jH] very difficult almost impossible. However, when the fluctuations 
of W are Gaussian we propose here an alternative method to compute AF, which can be useful 
even in cases where the JE and the CR could not hold. 

The paper is organized as follows. In section 2 we recall the JE and CR and we discuss the 
Gaussian case. In section 3 we describe the experimental setup and in section 4 the experimental 
results. In section 5 we show that for a Langevin equation the JE gives the exact value of the free 
energy difference for any path. In section 6 we discuss the limits of applicability of JE and CR 
and we conclude. 

2 The Jarzynski equality and the Crooks relation 

In 1997 pp Jarzynski derived an equality which relates the free energy difference of a system in 
contact with a heat reservoir to the pdf of the work performed on the system to drive it from A to 
B along any path 7 in the system parameter space. 

2.1 The Jarzynski equality 

Specifically, when A is varied from time t = to t — t s , Jarzynski defines for one realization of the 
"switching process" from A to B the work performed on the system as 



where z denotes the phase-space point of the system and H\ its A-parametrized Hamiltonian (see 
also 53] and section 2.4). One can consider an ensemble of realizations of this "switching process" 




(1) 
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with initial conditions all starting in the same initial equilibrium state. Then W may be computed 
for each trajectory in the ensemble. The JE states that 

AF = ~bx(expl-(3W]), (2) 

where (•) denotes the ensemble average, = fcsT with fcs the Boltzmann constant and T the 
temperature. In other words (exp [— /3Wdiss]) = 1, since we can always write W = AF + Wdi SS 
where Wdiss is the dissipated work. Thus it is easy to see that there must exist some microscopic 
trajectories such that Wdiss < 0. Moreover, the inequality (exp x) > exp (x) allows us to recover 
the second principle, namely (Wdiss) > 0> i- e - (W) > AF. From an experimental point of view the 
JE is quite useful because there is no restriction on the choice of the path 7 and it overcomes the 
above mentioned experimental difficulties. 

As we have already mentioned in the introduction, many proofs of the JE have been done. 
The simplest way to understand JE is to consider the perturbation theory JHlEj- This kind 
of approach has been used by Landau J3j but he considered only Gaussian distributions and 
he stopped the development to the second order. This technique can be generalized to get JE. 
Therefore we follow ref.113^ and we write the energy Hb(jp, q) of the system in contact with a heath 
bath at temperature T as: 

H B (p,q) = H A (p,q) + V(p,q,\(t s )) (3) 

where V is the energy introduced by the external driving A such that V(p, q, A(0)) = 0. By 
definition the AF is given by 

F B =F A -I3- I \nj PA exp {-(3V]dT (4) 

where p A (p,q) = exp [(3(Fa — HA{p,q))] is the equilibrium Gibbs distribution in A, that is the 
unperturbed state. Therefore 

F B = F A - I3- 1 In < exp [-f3V] > . (5) 
As = ?M- and H depends explicitly on time only by mean of A then 

Comparing eqElto eqQ]we see that V = W and eql^lis equivalent to eqEl 
2.2 The Crooks relation 

In our experiment we can also check the CR which is somehow related to the JE and which gives 
useful and complementary information on the dissipated work. Crooks considers the forward work 
Wf to drive the system from A to B and the backward work Wb to drive it from B to A. If the 
work pdfs during the forward and backward processes are Pf(W) and Pb(W), one has [21 E] 

= ex P W W - AF D = ex P \PW*»\- (7) 
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A simple calculation from Eq. Q leads to Eq. (JJJ. However, from an experimental point of view 
this relation is extremely useful because one immediately sees that the crossing point of the two 
pdfs, that is the point where Pf(W) = Pb( — W), is precisely AF. Thus one has another mean to 
check the computed free energy by looking at the pdfs crossing point W x . 



2.3 The Gaussian case 

Let us examine in some detail the Gaussian case, that is P(W) oc exp (- ■ J . In this case 



the JE leads to 



Qa 



AF = {W) L_jv ( (g) 



i.e. (Wdi ss ) = -y 51 > 0. It is interesting to notice that Landau J3j derived eqOO which can also 
be viewed as a consequence of the linear response theory as it has been shown by Hermans [Tfi] . 
It is easy to see from Eq. Q that if Pf(H^) and Pb(— W) are Gaussian, then 

AF= M*-(w)* t (9) 

and 

f3a 2 w = (W) t + (W) b = 2 (W diss ). (10) 

Thus in the case of Gaussian statistics AF and Wdiss can be computed by using just the mean 
values and the variance of the work W. 



2.4 The classical work and the AF computed by the JE 

Before describing the experiment, we want to discuss several important points. The first is the 
definition of the work given in Eq. JQ, which is not the classical one. Let us consider, for example, 
that A is a mechanical torque M applied to a mechanical system S, and —dH\/d\ the associated 
angular displacement 9. Then, from Eq. (QJ, one has 



W 



M6dt 



M6 



W c 



(11) 



where 



W c 



M6dt 



(12) 



is the classical work (we define the classical work with minus sign to respect the standard convention 
of thermodynamics) . Thus W and W cl are related but they are not exactly the same and we will 
show that this makes an important difference in the fluctuations of these two quantities. This 
difference between the W and W cl , has been already pointed out in ref.[T7|. 

The second important point that we want to discuss concerns the A_F computed by the JE 
in the case of a driven system, composed by 5 plus the external driving. The total free energy 
difference is 



AF = AFo - 



MO 



= AF - 



(13) 
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current J 




Figure 1: Schematic drawing of the oscillator 



where AFq is the free energy of S and <I> 



MB 



the energy difference of the forcing. The JE 



computes the AF of the driven system and not that of the system alone which is AF . This is 
an important observation in view of all applications where an external parameter is added to 3 
in order to measure AFq 0. Finally we point out that, in an isothermal process, AFq can be 
easily computed, without using the JE and the CR, if W is Gaussian distributed with variance 
a w cl - make the reasonable assumption that forward and backward work variances are equal.) 
Indeed the crossing point of the two Gaussian pdfs Pf(W cl ) and Pb(— W cl ) is 



(14) 



which by definition is just — AFo, i.e. — —AFq. Furthermore the dissipated work can be 

obtained from 

(Wdi SS ) = , (15) 



by definition. It should be noted that the equality 2 (Wdiss) = does not hold in the case of 

the classical work. 



3 Experimental setup 

3.1 The torsion pendulum 

To study the JE and the CR we measure the out-of-equilibrium fluctuations of a macroscopic 
mechanical torsion pendulum made of a brass wire, whose damping is given either by the vis- 
coelasticity of the torsion wire or by the viscosity of a surrounding fluid. This system is enclosed 
in a cell which can be filled with a viscous fluid, which acts as a heat bath. A brass wire of 
length 10 mm, width 0.75 mm, thickness 50 /im, mass 5.91 x 10~ 3 g, is clamped at both ends, 
hence its elastic torsional stiffness is C — 7.50 x 10 -4 Nmrad -1 . A small mirror of effective mass 

4.02 x 10~ 2 g, length 2.25 mm, width b\ = 7 mm, thickness a\ = 1.04 mm, is glued in the middle 
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(i) (ii) (iii) 




Figure 2: FDT check in oil. i) Real part 1/x. ii) Imaginary part of 1/x. iii) Noise spectrum 

the circles (o) correspond to the direct measurement of the noise and the dashed curve is \J (\0\ ) 
computed by inserting the measured response function \ m the FDT, Eq. I|18H 



of the wire, see Fig.^Ji), so that the moment of inertia of the wire plus the mirror in vacuum 
is / = 1.79 x 1CP 10 kgm 2 (whose main contribution comes from the mirror). Thus the resonant 
frequency of the pendulum in vacuum is /o = 326.25 Hz. When the cell is filled with a viscous 
fluid, the total moment of inertia is I c g = I + /fluid, where /fluid is the extra moment of inertia 
given by the fluid displaced by the mirror |Tj|. Specifically, for the oil used in the experiment 
(which is a mineral oil of optical index n = 1.65, viscosity 121.3 mPas and density p = 0.9p wa t e r 
at T — 21.3 °C) the resonant frequency becomes /o = 213 Hz. To apply an external torque M 
to the torsion pendulum, a small electric coil connected to the brass wire is glued in the back of 
the mirror. Two fixed magnets on the cell facing each other with opposite poles generate a static 
magnetic field. We apply a torque by varying a very small current J flowing through the electric 
coil, hence M ~ Amp J, where Amp is an amplification factor which depends on the size of the 
coil and of the distance of the magnets. This factor, which can be measured independently (see 
ref.[2l], is the largest source of error of our measurement. It is known with 3% of accuracy. The 
measurement of the angular displacement of the mirror 9 is done using a Nomarski interferometer 
|19H2()| whose noise is about 6.25 x 10~ 12 rad/VHz, which is two orders of magnitude smaller than 
the oscillator thermal fluctuations. An optical window lets the laser beams to go inside and outside 
cell. Much care has been taken in order to isolate the apparatus from the external mechanical and 
acoustic noise, see |22. for details. 

3.2 Equation of motion 

The motion of the torsion pendulum can be assimilated to that of a driven harmonic oscillator 
damped by the viscoelasticity of the torsion wire and the viscosity of the surrounding fluid, whose 
motion equation reads in the temporal domain 

Icsd+ f G(t-t')6{t')dt' + C6 = M, (16) 

J — oo 

where G is the memory kernel, which in the simplest case of a viscous damping is G(t ~ t') — 
vbit — t'). In Fourier space (in the frequency range of our interest) this equation takes the simple 
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form 



{-I cB lu 2 + C]9 = M, (17) 

where : denotes the Fourier transform and C = C + i[C +ojv] is the complex frequency-dependent 
elastic stiffness of the system. C and v are the viscoelastic and viscous components of the damping 
term. The response function of the system \ = 9/M can be measured by applying a torque with a 
white spectrum. When M = 0, the amplitude of the thermal vibrations of the oscillator is related 
to its response function via the fluctuation-dissipation theorem (FDT) ^2]. Therefore, the thermal 
fluctuation power spectral density (psd) of the torsion pendulum reads for positive frequencies 

(|^| 2 ) = ^lm*= ^ Z+P- . (18) 



3.3 Calibration and accuracy of the measurements 

To calibrate the system we check whether FDT is satisfied by the measurements. We apply a 
current J with a white spectrum. We measure the transfer function M/9 = 1/x = [— 7 e ff u> 2 + C + 
i(C" +U3v)]. As an example we plot in figEl) the measured real part of 1/x and in ii) the imaginary 
part. The parabolic fit of the real part and the polynomial fit of the imaginary part (v may depend 
on frequency) allow us to estimate all the parameters of the oscillator. The main source of error 
is given by the amplitude of the applied torque which is known with about 3% accuracy. Once \ 
is known the power spectrum of 9 with M = can be computed from eqEl We plot in Fig.^iii) 
the measured thermal square root psd of the oscillator. The measured noise spectrum [circles in 
Fig.|3{iii)] is compared with the one estimated [dotted line in Fig-EJiii)] by inserting the measured 
X in the FDT, Eq. 1I18J1 , The two measurements are in perfect agreement and obviously the FDT 
is fully satisfied because the system is at equilibrium in the state A where M — (see below). 
Although this result is expected, this test is very useful to show that the experimental apparatus 
can measure with a good accuracy and resolution the thermal noise of the macroscopic pendulum. 



4 Experimental results 



Now we drive the oscillator out of equilibrium between two states A (where M = 0) and B 
(where M = M m = const ^ 0). The path 7 may be changed by modifying the time evolution 
of M between A and B. We have chosen either linear ramps with different rising times r, see 
Fig-Eli); or half-sinusoids with half-period r. In the specific case of our harmonic oscillator, 
as the temperature is the same in states A and B, the free energy difference of the oscillator 



alone is AF = AU 



\C9 2 



2C 



AFn 



M 



i.e. for an harmonic 



whereas AF 

potential AF — —AFq. Let us first consider the situation where the cell is filled with oil. The 
oscillator's relaxation time is given by the inverse of the line width of the equilibrium fluctuation 
spectrum, see Fig.mii), that is r re i a x = 23.5 ms. We apply a torque which is a sequence of linear 
increasing /decreasing ramps and plateaux, as represented in Fig.^i). We chose different values of 
the amplitude of the torque M [22.1 11.9, 6.1, 4.2 and 1.2 pNm] and of the rising time r [199.5, 
20.2, 65.6, 99.6, 2.5 ms], as indicated in Tabled [cases a). . .g)] (for the case f) and g) also C has 
been changed). Thus we can probe either the reversible (or quasi-static) paths (r 3> r ro i ax ) or the 
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Figure 3: Oscillator immersed in oil [case (a)]: (i) Applied external torque, (ii) Induced angular 
displacement, (iii) its psd, (iv) its pdf, (v) Injected power computed from the Jarzynski definition 
W = —M9, (vi) Injected power computed from the standard definition W cl = —M9 



irreversible ones (r -C r re iax)- We tune the duration of the plateaux (which is at least 4r rc iax) 
so that the system always reaches equilibrium in the middle of each of them, which defines the 
equilibrium states A and B. We see in Fig.EJii), where the angular displacement 9 is plotted as a 
function of time [case a)], that the response of the oscillator to the applied torque is comparable 
to the thermal noise spectrum. The psd of 9 is shown in Fig.OJiii). Comparing this measure with 
the FDT prediction obtained in Fig.^Jii), one observes that the driver does not affect the thermal 
noise spectrum which remains equal to the equilibrium one. Moreover we plot in Fig.OJiv) the pdf 
of the driven displacement 9 shown on Fig.OJii), which is, roughly speaking, the superposition of 
two Gaussian pdfs. From the measure of M and 9, the power injected into the system W can be 
computed from the definition given in Eq. Q, that in this case is W = —M0. Its time evolution, 
shown in Fig.^Jv), is quite different from that of the classical power W cl = —M0, whose time 
evolution is plotted in Fig.|^vi): W is non-zero only for M/0 and vice- versa W cl ^ only for 
M ^ 0. From the time series of W we can compute from Eq. the forward and the backward 
works, Wf and Wb, corresponding to the paths A — » B and B — > A, respectively. We also do the 
same for the classical work. We then compute their respective pdfs Pf(W / ) and Pb(— W). These 
are plotted on Figs.2Ji,iv) where the bullets are the experimental data and the continuous lines 
their fitted Gaussian pdfs. In Fig.Q] the pdfs of W and W cl cross in the case a) at (3W ~ —23.5, 
and in the case c) (3W ~ —6.1. These values correspond to AF = —^f = ~ Ai 7 !)- We find that 
this result is true independently of the ratio r/r r eiax and of the maximum amplitude of |M|, M m . 
This has been checked at the largest M m and the shortest rising time t allowed by our apparatus. 
Indeed torques with amplitudes larger than 25 pN m and rising time shorter than 2 ms introduce 
mechanical noises which can be higher than thermal fluctuations and the check of the JE becomes 
impossible (see discussion in the next section). The measurement at very large M m and very 
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Figure 4: Case a): (i) P f (W) and P b (-W), (ii) P f (W cl ) and P b (-VK cl ); Case c): (iii) P f (W) and 
Pb(— W)j (iv) Pf(W^ cl ) and Pb(— W cl ) (experimental forward and backward pdfs are represented 
by o and □ respectively, whereas the continuous lines are Gaussian fits) 



short t is shown in Fig.0 Also in this case we see that the shape of the thermal noise spectrum 
is not perturbed by the driving. The pdfs remain Gaussian but the distance between (W) { and 
(W) b is larger and the relative variance much smaller than at low amplitude. The crossing point 
of the pdfs of W occurs at a value where the statistics is very poor. However the two Gaussian 
fits crosses at AF. The pdfs of W cX also crosses at the right values. The experimental results 
are summarized in Table [H where the computed AU = -^f- is in good agreement with the values 
obtained by the crossing points of the forward and backward pdfs, that is AF X + $ for P(W) and 
—AW^ for P(W /cl ). Finally inserting the values of Wf and Wb in Eq. ® we directly compute 
AFf and AFb from the JE. As it can be seen in Table the values of AFo obtained from the 
JE, that is either —(AFf + <£>) or — (AFt, + <&), agree (see sec. 3. 3) with the computed AU within 
experimental errors, that are about 5% on AU (see sec. 3. 3). The JE works well either when 
t ^> T rc i ax or in the critical case f) and g) where t <C r ro i ax . The other case we have studied 
is a very pathological one. Specifically, the oscillator is in vacuum and has a resonant frequency 
fo — 353 Hz and a relaxation time r re i a x = 666.7 ms. We applied a sinusoidal torque whose 
amplitude is either 5.9 x 10~ 12 or 9.4 x 10 -12 Nm [cases h) and i) in TableQ] respectively]. Half 
a period of the sinusoid is r = 49.5 ms, much smaller than the relaxation time, so that we never 
let the system equilibrate. However, we define the states A and B as the maxima and minima 
of the driver. Surprisingly, despite of the pathological definition of the equilibrium states A and 
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Figure 5: (i) Time evolution of 9 for the case g) in Table (ii) Its psd y {\9\ ), (iii) Pf(W / ) and 
Pb(— W), (iv) Pf(W cl ) and Pt>(— W ci ) (experimental forward and backward pdfs are represented 
by o and □ respectively, whereas the continuous lines are Gaussian fits) 



B, the pdfs are Gaussian and the JE is satisfied as indicated in Table Moreover, this happens 
independently of M m and of the critical value of the ratio r/r rc i ax -C 1. Finally, we indicated in 
Table the value AFq which is the free energy computed from the JE if one considers the "loop 
process" from A to A (the same can be done from B to B and the results are quantitatively the 
same). In principle this value should be zero, but in fact it is not since we have about 3% error in 
the calibration of the torque M and on C. 

5 Jarzynski equality and the Langevin equation 

In the previous section we have seen that even at very large driving torque, two properties of 
the systems remain unchanged, specifically: the thermal noise amplitude and statistics are not 
modified by the presence of the large forcing, that is FDT is still valid and the fluctuation pdf 
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Table 1: Free energies of cases a). . . g) defined in the text (the values of M m are in pNm). AU = 

^f- is the computed expected value. Notice that C — 7.5 x 10~ 4 Nmrad -1 for cases a)-e), 
C = 5.5 x 10~ 4 Nmrad -1 for cases f)-g) and C = 1.6 x 10~ 4 Nmrad -1 for cases h)-i) 



remains Gaussian. As we have already discussed, the fact that FDT is still valid is easily understood 
by comparing the noise spectrum of figGJh) without driving with those with driving in figE^ii) and 
figlSU). We clearly see that the driving does not change the shape of the noise spectrum, once the 
pics corresponding to the driving have been subtracted. It is also easy to show that the pdf Pd(0) 
of 9, with the driving, is the convolution product of the pdf P{&), without driving, times P(6) 
which is the mean response of to the driving torque. As an example we show in figEK) P{&) 
without driving. The PDF P{9) of the mean response of 9, measured for the case g) of Table is 
plotted in figEb). This PDF is close to the sum of two delta functions. In figEfc) we show, for the 
case g) of Table m the directly measured Pd{9) (circles) and that computed by the convolution of 
P{9)) with P{9). The agreement is excellent. This result shows that the statistical properties of 
the thermal noise are those of equilibrium even when the system is driven very fast. 
Starting from these experimental evidences, we can show that for the Langevin equation the JE 
is always satisfied independently of M m and r. Let us consider the equation for the harmonic 
oscillator, Eq. iflfijl . which well describes our experimental system. In the case of a viscous damping 
we rewrite Eq. Il l ( ill as 

I eS 9 + v6 + C9 = M + r), (19) 

where r\ is the thermal noise amplitude. For M(t) we consider the kind of waveform used in the 
experiments: 

M(t) = for 0<t<r, (20) 

T 

= M m for t > t. (21) 
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Figure 6: (i) PDF of 9 when the driving is off. (ii) PDF P(9) of the mean response of 9 to the 
torque, measured for the case g) of Table (iii) PDF Pd{9) of 9 measured for the case g) of Table 
Ql(o) and computed by the convolutions of the PDF plotted in (i) and (ii) (continuous line). 



5.1 The overdamped case 

In order to compute AF using the JE when the system is driven from A (M = 0) to B (M = M m ^ 
0), we first consider the overdamped case when the inertial term is negligible, that is I e g — 0: 

9 + T^ 1 9 = (CT y 1 (M + r ] ), (22) 

where r = vjC. If 77 is a thermal noise, then when M = the spectrum of the thermal fluctuations 
89 of 6 can be computed from FDT: 



~,2 N 4fc B Tr 



As a consequence the autocorrelation function of 9 on a time interval 5t is 



(23) 



W^e^f-M). (24) 



C 1 V r 

The pdf of 9 is Gaussian. The work to drive the system from A to B computed using Eq. 
becomes in this case 

w ^_Mrn r edt (25) 

T Jo 

Thus to compute W we need only the solution of Eq. 112211 for < t < t. If we neglect the noise, 
then the mean solution is 



tC 



t + T CXp ( ) - T 

T~0 



for < t < t. (26) 



We now consider that in the experiment the statistical properties of the thermal fluctuations 59 
are not modified by the driving. Therefore 9 can be decomposed in the sum of an average part 9 
given by Eq. 12611 plus the fluctuating part 59, that is 

9 = 9 + 59. (27) 
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As a consequence the work can also be decomposed in a similar manner 

W = W + SW 



T 







[ 6dt + 


f 69 dt 


Jo 


Jo 



(28) 



As the integral of a Gaussian variable is still Gaussian then the fluctuations of W remain Gaussian 
too. As a consequence, to compute AF we can use Eq. where (W) is straightforward computed 
using Eqs. ® and QBJ 



/W\ W Mm [ T HA+ 1 fM m Yl (T-T ) 2 2 / T\ TO 2 

We now have to compute erjv = (($W) 2 ) = (M m /r) 2 (y 2 (r)) where 



>dt. 



The variance of y can be computed j2l] taking into account that 

<2/ 2 W> = / / flfl(ti-t 2 )dtidta. 
Jo Jo 

Using this equation and Eq. ll23jl we get 



2 



-Mir 
r 



(y 2 W> = 



2k B TT [M, 
C~ 



t - t + t exp - 



T() 



(29) 



(30) 



(31) 



(32) 



Taking into account that fluctuations of W are Gaussian, we replace the results of Eqs. l(29"jl and 
il3*2")l in Eq. © , and finally we get 



AF = (W) 



4 
2k B T 



M 2 
2C ' 



(33) 



that is the expected value. It is important to notice that this equation gives the exact result 
independently of the rising time of the external applied torque. This is important because it 
shows that the in cases where the conditions a) and b) are verified, the JE gives the right result 
independently of the path to go from A to B, which can be a very irreversible one. 

5.2 The harmonic oscillator 

We may now repeat the calculation for the harmonic oscillator of Eq. iflQll . In such a case let us 
introduce the following notations 



and 



cos ip = 



2-Zeff 



\/a 2 + tp 2 



a 2 + t/j 2 



sin ip 



C_ 



-J— 

y/a 2 + -0 2 



(34) 
(35) 
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With this notation the mean solution of Eq. ljl9H in absence of noise, with initial conditions 9(0) = 
0(0) = 0, is 



Mm 
Tip 



exp (—at) sin (^)i + 2<£>) + ip t — sin 2tp 



for 



< t < r. 



The correlation function of the thermal fluctuations 89 become in this case 



Re(6r) = 



C siiup 



exp (— a|5r|) sin (?/>|5t| + f). 



(36) 



(37) 



We now proceed as in the overdamped case and we compute (W) and a^. Using Eq. ll-'ilit and 
Eq. ffify we get 



siny 



exp (— ar) sin (ipT + 3ip) — sin3</? — 



(V>t) 2 

2sin< / 9 



2^r cos (/? 



(38) 



To compute the variance of W, we insert Eq. lTT7jl in Eq. (131)) and we obtain 



2 



M„ 

T 



(y 2 (r)) 



M m \ 2ksT sin ip 

ex P \~ aT ) sm (v 7 " + 3y>) — sm 3o? + 2f/'T cos 93 . (39) 

T / CV' 

As the fluctuations of W are Gaussian, we insert Eqs. jSEt and (l39l) in Eq. JSJ, and we obtain 



ow Mi 



AF = (W) ^— = 

v 7 2fc B F 2C ' 



(40) 



which is the expected results. Notice that in this case too the result is independent on the path. 
Thus the JE gives the right result of AF for the Langevin equation with an harmonic potential. 



6 Discussion and Conclusions 
6.1 Experimental limits of JE and CR 

Before concluding we want to discuss the limits on the use of JE and CR in an experiment where 
the dynamics can be either exactly described or well approximated by a Langevin equation. From 
the values of (W) and of aw, computed in sect. 4 for the overdamped case and for the harmonic 
oscillator case, we see that if r is kept constant then (W) oc whereas aw oc M m . Furthermore 
as 2 AF = (W) f — (W) b is proportional to M^, this means that the distance between the maxima 
of Pf(W) and Pb(— W) increases with Mf n and the relative width |ervp/(W)| of P(T / F) decreases as 
X/M m . As a consequence, the probability of finding experimental values of W close to the crossing 
points also decreases as 1/M m . This effect has been seen on the experimental P(W) plotted in 
Fig.0 From a practical point of view, this means that when AF >• fc^T the pdfs will never cross, 
for reasonable values of the work pdf. As a consequence the experimental test of JE and CR in 
systems with a strongly irreversible driving force, as proposed in ref.^T], will be very difficult to 
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test. Indeed at such a large and fast driving, the JE and CR may fail, either for deep theoretical 
reasons or simply because the statistics is too poor. At the moment we are not able to answer to 
this question. Furthermore, for real (macroscopic) systems it is reasonable to think that when the 
external noise becomes much larger than the thermal noise, the JE cannot be used. However, if 
the pdfs of W cl remain Gaussian, then the crossing point gives the right result. 

6.2 Summary and discussion of the results 

In conclusion we have used a driven torsion pendulum to test experimentally the accuracy of 
JE and CR. By varying the amplitude and the rising time of the driving torque of about one 
order of magnitude, we clearly demonstrate the validity and the robustness of the JE and CR in 
an isothermal process, at least when the work fluctuations are Gaussian and when the harmonic 
approximation is relevant for the system. We have checked the generality of the results on a driven 
Langevin equation which well describes the dynamics of the oscillator. Using the experimental 
observations that the equilibrium properties of the thermal noise are not modified by the driving 
force and that the force fluctuations are Gaussian, we have shown that the JE gives the right 
result independently of M m and r. Unfortunately these results do not fully alight the theoretical 
debate, because our conditions are close to those pointed out in Ref. f° r the validity of the JE. 
Recently, Ritort and coworkers have used the JE and CR to estimate AF in an experiment of RNA 
stretching where the oscillator's coupling is non-linear and the work fluctuations are non-Gaussian 
|25j . It would be interesting to check these results on a more simple and controlled system. We 
are currently working on the experimental realization of such a non-linear coupling, for which 
AF 7^ — AF . 

We have also shown both analytically and experimentally that |ervv/(WOI decreases as 1/M m . 
This observation makes the practical use of the JE and CR rather unrealistic for very large M m 
as the statistics needed to get a reliable result will be very large. This means, as we have already 
discussed at the end of the previous section, that the JE cannot be applied for macroscopic systems 
when the external noise becomes much larger than the thermal noise. 

Going back to the estimate of AF we have seen that the more accurate and reliable AF 
estimator is given by the crossing points AF X and W 1 ^, because they are less sensitive to the 
extreme fluctuations which may perturb the convergence of the JE. Starting from this observation, 
we propose a new method to compute AF in the case of Gaussian fluctuations of W cl . Indeed, for 
Gaussian pdfs W$ remains an excellent estimator even in cases where the JE and the CR could 
not hold, for example when the environmental noise cannot be neglected. 

Finally we want to stress that our results, although limited to the Gaussian case, show that it 
is possible to measure tiny work fluctuations in a macroscopic system. As a consequence it opens 
a lot of perspective to use the JE, the CR and the recent theorems on dissipated work (see for 
example 26j) to characterize the slow relaxation towards equilibrium in more complex systems, 
for example aging materials such as glasses or gels [27 . 

The authors thank L. Bellon, E.G.D. Cohen, N. Gamier, C. Jarzynski, F. Ritort and L. Rondoni 
for useful discussions, and acknowledge P. Metz, M. Moulin, F. Vittoz, C. Lemonias and P.-E. Roche 
for technical support. This work has been partially supported by the Dyglagemem contract of 
EEC. 
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